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ABSTRACT 

We perform a strong- lensing analysis of the merging galaxy cluster MACS J0416. 1-2403 (M0416; 
z = 0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 
23 background sources in the range 0.7 < z p hot ^ 6.14 including two probable high-redshift dropouts, 
revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~ 100" (z s ~ 2). Compared 
to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area 
is not particularly large (~ 0.6 □'; z s ~ 2), the number of multiple images, per critical area, is 
anomalously high. We calculate that the observed elongation boosts the number of multiple images, 
per critical area, by a factor of ~ 2.5 x, due to the increased ratio of the caustic area relative to 
the critical area. Additionally, we find that the observed separation between the two main mass 
components enlarges the critical area by a factor of ~ 2. These geometrical effects can account for 
the high number (density) of multiple images observed. We find in numerical simulations, that only 
~ 4% of the clusters (with M V i r > 6 x 10 14 /i _1 M Q ) exhibit as elongated critical curves as M0416. 

Subject headings: dark matter, galaxies: clusters: individuals: MACS J0416. 1-2403, galaxies: clusters: 
general, galaxies: high-redshift, gravitational lensing: strong 



1. INTRODUCTION 

Due to their high, projected surface mass densities, 
galaxy clusters magnify and distort background objects, 
forming natural gravitational lenses in the sky. The lens- 
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ing and magnification effects generally increase towards 
the central region of the cluster, where the projected 
mass density is often high enough to form multiple im- 
ages of the same background source, depending also on 
the angular diameter distances involved (for reviews see 
Bartelmann||2010MKneib k Natara"jan"||2011 1. 

The lensing efficiency of galaxy clusters (e.g. the num- 
ber of multiple-images generated) is also related to other 
factors, such as the ellipticity, amount of substructure 
and its distance fro m the center, and degree of relax- 
ation or merger (e.g. Meneghetti et al. 2003[ ). For ex- 



ample, the critic al area grows with the co ncentration of 
the cluster (e.g. Sadeh k Rephaeli 2008), so it is clear 
that massive and more concentrated clusters should show 
more multiple images. The "over-concentration" prob- 
lem in which lensing-selected clusters are found t o have 
high concentrations (and large Einstein radii, e.g . Com- 



erford k Natarajan|20l^|Broadhurst et al.l2008||Broad 

hurst k Barkana 2008), is often attributed to a lensing 



selection bias towards higher concentrations of triaxial 
clus ters preferentially aligned with the line of sight (see 
also|Hennawi et aL 2007; Sc reno et al.|2010| |Oguri et al. 
2012a[|Okabe et al.||2010 and references therein) . 

On the other hand, recent efforts show that there ex- 
ists a s econd class of prominent strong lenses. Sereno fc 
(2012) showed that in a tria xial lensing analysis 



Zit 



the 12 MACS clusters at z > 0.5 (|Ebeling et al 20071 



* Mass models are publicly available at: ftp://wise- 
ftp.tau.ac.il/pub/adiz/M0416/ 



have relatively low concentrations, despite the many mul- 
tiple images uncovered in their fields. They suggested 
that since most of these clusters are not yet relaxed, the 
amount of substructure in their centers as well as their 
higher redshift than most previously known lenses, turn 
them into highly-magnifying lenses. The critical curves 
of the several subclumps are merged together into a big- 



2 



Zitrin et al. 



ger lens, whose overall mass profile in the central part is 
often shallower thus boosting the magnification. 

Indeed, various clusters around z ~ 0.5 show promi- 
nent lensing features (giant arcs and many multiple im- 



20 clusters are X-ray selected to be dynamically relaxed). 
We found no record of a published strong l ensing (SL) 



ages; IZitrin et al.||2012a|bp. Interestingly, the largest 



gravitationa l lens, MACS J0717.5+3745 (z = 0.55; Zitrin 



et al.|200"9a ), is formed by s everal merging clumps possi 



re giant arc 
1.896, which 



bly at the tip of a filament (Limousin et al. | 2012 Jauzac 



et al. 



2012) 



showed that 
0.54) is an excellent "mag- 



Zitrin & Broadhurst 
MACS J 1 149.5+2223 (z 
nifying glass" in the sky due to its shallow inner mass 
profile. T his cluster is now known to magnify a galaxy 
at z ~ 10 flZhenget al.|2012[ ). Recently, [Coe et al.|p012 1 
uncovered the highest redshift galaxy known to date 
at z ~ 10.8, multiply-imaged by MACS J0647.7+7015 
(z = 0.59), a nother complex cluster expected to be h ighly 



magnifying ([Zitrin et al.||2011 [ |Postman et al. 



2012). 



|Meneghetti et al.| ( |2U03p found that the cluster lensing 
cross section for giant arcs grows rapidly with elliptic- 
ity. A typical ellipticity of e = 1 — b/a = 0.3 entails an 
order of magnitude increas e in the lensing cross section 



(see also Torri et al. 
increasing ellipticity 

and enclose a growing area. |Meneghetti et al.| (|2003[) 



2004J), for example, because with 
he caust ics stretch, develop cusps 



also found, that ap proaching subclumps boos t the lens 
ing cross sections. Meneghetti et al. (2007) examined 
the arc sensitivity to cluster ellipticity, asymmetries, and 
substructures, and found that these contribute, respec- 
tively, ~ 40%, ~ 10%, and ~ 30%, t o the lensing cross 
section. Recently, |Redlich et al.| ( |2012| ) found that cluster 
mergers can enhance the lensing cross section by typically 
~ 30 — 50%. It is therefore expected that merging, sub- 
structured and elongated clusters, should also constitute 
prominent strong lenses. 
Here, we present the lensing analysis (§2) of the merg- 



analysis of M0416, but Christensen et al. (|2012|) pub- 
lished a measured spectroscopic redshitt for t 
north of the BCG (systems 1 & 2 here) 
we use in our analysis. 

Throughout we adopt a concordance ACDM cosmology 
with (Q m0 = 0.3, O\ = 0.7, h = 0.7), where 1" = 
5.53 kpc at the adopted redshift for M0416, z = 0.42 
( |Christensen et al.||2012[ ). 

2. STRONG LENSING ANALYSIS OF M0416 

We use two complementary modeling techniques to 
construct mass models for M0416. These are then com- 
pared to each other and to our new weak lensing (WL) 
analysis. 

Method 1 : The first method we use follows the pre- 
scription of Zitrin et al. (2009b), with several modifica- 
tions recently implemented tor speed and a wider choice 
of priors. First, instead of power-law profiles tr adition- 
ally used to represent the galaxies in ou r method ( Broad- 
hurst et al. |2005| Zitrin et al. 2009b), we model each 
galaxy with a P seudo-Is othermal Elli ptica l Mass Distri- 
bution (PIEMD; |Kassiola fc Kovner|1993| ), adopting the 
exact formulation trom |Jullo et al.| ( f2007[ ): 



^(^) 1/4 



(1) 



r =r * 

' core ' core\ L* ) 
Tcut — r ctI t(jT") i 

where r core is the core radius, r cut the cut-off radius, and 
do the velocity dispersion. The total mass of a subhalo 
then scales as: 



M = (ir/G)(** ) 2 r* cut (L/L*) 1/2+a , 



(2) 



ing cluster MACS J0416. 1-2403 (hereafter M0416; Mann where L* is the typical luminosity of a galaxy at the 



& Ebeling |2012 ), performed on recent Hubble Space 
Telescope (HST) imaging in 16 bands from the UV to 
the near-IR to a total depth of ~ 20 orbits, as part of 
the Cluster Lensing And S upernova survey with Hub- 



cluster redshift, and r* 



parameters (Julio et al. 2007) 



and cr? are its PIEMD 



The PIEMD representati on is used in many model 



ble (CLASH) program (see |Postman et al.||2012| ). The 
CLASH pipeline uses the 16-band observations to de- 
rive p hotometric redshifts for each arc, via th e BPZ pro- 



ing methods (e.g. |Halkola et al.||2006 




Julio et al. 


2007 


|Richard et al.| 2010| |Uguri et al.||2012 


b), and eases the 



gram (Bemtez et al. 



2004 Coe et al. 2006), used here 



to constrain the model. As we shall see, the cluster ex- 
hibits high elongation, in part as a result of the possible 
merger. We test, by a simple semi-analytical simulation 
(§3), the expected increase in lensing efficiency with el- 
lipticity and merger stage, to see if these may account 
for the many multiple images in M0416, per its critical 
area. For this work, we define the lensing efficiency sim- 
ply as the number of multiple images per critical area, 
which we also refer to as the number density of multiple 
images. This ratio, naturally, scales with the ratio of the 
caustic area and the critical area. 

M0416 w as listed as a MACS cluster (MAssive Cluster 
Surve y, see Ebeling ct al. 20l0l) due to its X-ray bright- 
ness. |Mann fe Ebeling |2012| classified it as merging 
based on the double-peaked X-ray structure, where the 
southern peak is offset by a few arcseconds from the cor- 
responding (second) BCG, as expecte d in mergers (e.g. 
Bradac et al.|2006 Merten et al.|20lT ). Based also on its 
predicted Einstein radius, M0416 was designated as one 
of the five "high-magnification" CLASH clusters (where 



a second change, this mass distribution is now smoothed 
with an elliptical Gaussian (eGaussian), instead of our 
traditional smoothing by a 2D polynomial spline, and 
therefore eliminates the need for an external shear as 
the ellipticity is now directly introduced into the mass 
distribution. The superposed galaxies' contribution and 
its eGaussian-smoothed map - are then added with a 
relative weight that is a free parameter. This method 



therefore includes 6 free basic parameters: 



for 



the PIEMD galaxy models (eqs. [T]and[2]); o~mja and 
&mna, the widths of the eGaussiaiiKernelalong the ma- 
jor and minor axes; <j>, the angle of the major axis in the 
eGaussian kernel; and the relative fraction of the galaxies 
component from the total mass (the remaining fraction 
is contributed by the smooth DM component). 

Method 2: The second method we use adopts the fol- 
lowing parametrization. Galaxies are modeled each as 
PIEMD scaled by its light as above (Method 1; eq. 
[I). To represent the DM we implement an elliptical 
JNFW (eNFW) halo. Since M0416 is clearly a com- 
plex, likely merging system, the first few modeling at- 
tempts introduced the need to add a second eNFW 
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Fig. 1. — CLASH/HST 16- band color-composite image of M0416, with multiple images numbered. Lensed candidates are marked with 
"c", and "p" stands for predicted location. Overlaid in blue(white) is the critical curve for z s = 1.896, corresponding to the giant arc 
[system 1&2], from the eNFW(eGaussian) based model (see §2). In red(green) we plot the corresponding caustics. 



halo to the model. In total, then, the model includes 
the galaxy component described by the superposition 
of all PIEMD representations, and two eNFW halos 
(where the ellipticity, defined here throughout as e = 
[MJA - MNA)/(MJA + MNA), is directly introduced 
into the NF W mass distribution via the transformation 
r — > r e = y/[x/{l + e)] 2 + [y/(l — e)] 2 ). We maintain 
the eNFW halos centered on the first and second bright- 
est cluster members, respectively. This keeps the number 
of free parameters as low as possible and yields an excel- 
lent fit. Note that we tried to construct a model while 
allowing the (southern halo) center to vary but did not 
find a significantly better solution. This parametrization 
therefore consists of 10 free basic parameters: r* ut , cTq for 
the PIEMD galaxy models; the scale radius r s and the 
concentration parameter c V i r , as well as the ellipticity 



and its position angle, for each of the two eNFW halos. 

Using a preliminary model with the |Zitrin et al. 
( 2009b ) method, and the two iteratively-improved mod- 
els above, along with a complementary examination by 
eye, we iteratively matched together multiply-imaged 
systems. The best fit solution in both methods is ob- 
tained via a long (several dozens of thousands steps) 
Monte-Carlo Markov Chain (MCMC) minimization, us- 
ing several chains. We note that in some chains, includ- 
ing the final ones used here, the redshifts of some of the 
multiple-systems were left as free parameters (with flat 
priors) to be optimized. Also, in both final chains for 
methods 1 and 2, we left the three BCGs to be freely 
weighted and optimized by the MCMC. For the two main 
BCGs, we fix the ellipticity and pos ition angle to the 
parameters derived using SExtractor (Bertin & Arnouts 
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Fig. 2. — Projected mass density profile of M0416. The red curve 
shows the resulting profile and la errors from the two eNFW halos 
model; and the black curve shows the resulting profile from the 
eGaussian-smoothing model (jj2]for details). The two models are 
similar in the range where multiple- image constraints are available 
(< 1'), and highly consistent with independent Subaru WL analysis 
data (blue squares and errorbars). 



1996). Throughout, we fix for the PIEMD, r* core = 0.3 



kpc, use an L* value equivalent to an absolute magnitude 
of Mp 814W = —20.113, and assume a constant mass-to- 
light ratio {a = 0.5, eq. pj. The F814W-F475W color 
was used to extract the recFsequence, where we used here 
the 122 brightest (F814W) members. 



2.1. 



Results and comparison of the two models 



In total, we found 70 multiple images and candidates 
of 23 background sources (Table [I] Fig. [j). All images 
not marked as candidates were used as constraints for the 
model: 34 images of 13 sources. For seven of these the 
redshifts were left to be optimized by the models. Two 
of the multiply-imaged sources seem to be high redshift 
"dropouts": system 6 at z ~ 6.5, and candidate system 
20 at z ~ 5, although the redshift constraints on the 
latter are poor. We leave the detailed examination of 
these two objects for future work. 

Both models yield very similar critical curves, ex- 
cept for a small region of a few arcseconds discrepancy 
where there are only systems that were not used as con- 
straints (or their redshifts left free). In addition, the 
mass profiles (centered on the midpoint between the first 
and second brightest members, RAj 2 ooo=04:16:08.38, 
DECj2ooo=-24:04:20.80) of the two models are almost 
identical throughout the range where multiple images 
are observed (r < 1'), see Fig. [2j The two models are 
also in excellent agreement wit h independent, color-color 
selected background galaxies ( |Medezinski et aL 20101 
WL measurements from BR c z Subaru data (Fig, pj)' , 
using the Bayesian method of Umetsu et al. ( 2011a]b| ) 
that combines tangential-distortion and magnirication- 
bias measurements in a model-independent manner, ef- 
fectively breaking the mass-she et degeneracy. The W L 
analysis pipeline is described in Umetsu et al. (2012). 

The image-plane reproduction rms(x^) is 1.89" (56. 67) 
and 1.37" (29.67) for the eGaussian and eNFW mod- 
els, respectively. The rms increases slightly to 2.40"and 
2.11", respectively, when including all images and can- 



didates and not only those used as constraints. For 
the x 2 we used a positional error of a pos = 1.4", 
which was found to be a reasonable value accounting for 
larg e-scale structure and matter along the line-of-sight 
(see |Zitrin et al. 2012b| and references therein). The 
multiple-images input comprises 35 constraints, where 
the number of DOF is 26 and 22 for the eGaussian 
and eNFW models, respectively, yielding correspond- 
ingly x 2 /DOF=2A8 and % 2 /£)OF=1.35. Both models 
show a critical area (A) with an effective Einstein ra- 
dius of 6 e = ~ 26 ± 2" (z s = 1.896), enclosing 

M(< 9 e ) = 1.25 ± 0.09 x 10 14 M Q . Comparing the two 
methods, the eNFW model has a somewhat better fit to 
the data - and higher flexibility to match it. The eGaus- 
sian smoothed model fit is somewhat worse, following 
more rigourously the light-traces-mass assumption, but 
the fact that it physically matches multiple images a pri- 
ori (without needing to accurately constrain the fit first), 
is remarkable evidence for its credibility. 

The number density of multiple images and candidates 
uncovered, i.e. number per the given critical area (70 
images of 23 sources over ~ 0.6 □' for z s ~ 2), is a few 
times higher than other prominent known lensing clusters 
(for similar background redshifts). A1689 (z = 0.19) 
for example, one of the largest lenses known, shows 135 
images of 42 sou rces, over a critical area of ~ 1.8 □' 
( |Coe et al.||2010| ). A1703 (z = 0.28) has more than 50 
multiple im ages of 17 sources kn own, and a critical area 
of ~ 0.8 □' ( (Richard et al][2009| . MACS J 0717.5+3745 
(z = 0.55), the largest gravitational lens (Zitrin et al. 



2009a ) , has abou t 60 multiple images known trom 18 
sources (see also Limousin et al. 2012), and a critical 
area of ~ 2.64 □'. 

3. EXPECTED EFFECT OF ELLIPTICITY AND MERGER 

We now wish to test the effects of ellipticity and merger 
on the observed number of multiple images. 

We start by producing a fiducial eNFW model at 
z = 0.42, with point sources planted behind it [z s = 2), 
every 45 kpc on a grid. We then increase the elliptic- 
ity (e = (MJA - M NA) / (M J A + MNA)) from 0.0 to 
0.8, and count the number of multiple images generated. 
Note that throughout we always count all multiple im- 
ages formed. For each configuration we also measure the 
critical area for normalization, so the effect of ellipticity 
can be extracted per given critical area. 

In Fig. [3] we plot the resulting increase in the num- 
ber of multiple images with lens ellipticity, per critical 
area, for various combinations of NFW parameters. A 
clear correlation is observed, so larger numbers of multi- 
ple images are generated by higher ellipticities (but the 
amplitude may vary with masses and distances). To as- 
sess the effect in M0416 (for a fixed axis ratio of ~ 5.1), 
we choose different combinations of C200 and r s that yield 
comparable critical area to that of M0416. The observed 
elongation results in a ~ 2.5x increase in the number of 
multiple images compared with the circular case, for the 
given critical area. Therefore, the observed elongation 
explains why M0416 has a few times higher multiple- 
image number density than other typical lensing clusters 
(in comparable HST imaging). In reality, the increase 
in the number density of multiple images with elliptic- 
ity is of course, a more complex function and the exact 
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TABLE 1 

Multiple Images and Candidates found by our method. 



Arc ID 


ra(J20oo.o) 


DEC(J2000.0) 


phot-s [z min -z max ] 


Z NFW 


A NFW " 


'Gauss 


A Gauss " 


comments 


1.1 
1.2 
1.3 
2.1 
2.2 
2.3 


04 
04 
04 
04 
04 
04 


16:09.780 
16:10.435 
16:11.365 
16:09.884 
16:10.321 
16:11.394 


-24:03:41.73 
-24:03:48.75 
-24:04:07.21 
-24:03:42.77 
-24:03:46.93 
-24:04:07.86 


1.788 [1.541-1.879] 
2.482 [2.379-2.560] 
2.555 [2.442-2.675] 
1.788 [1.541-1.879] 
1.846 [1.796-1.998] 
1.928 [1.806-2.239] 


2.01 [1.66-2.69] 
2.00 [1.71-2.13] 


1.6 
0.9 
0.3 
1 

0.5 
0.7 


2.10 [1.67-2.17] 
2.24 [1.73-2.24] 


1.5 
2.3 
0.9 
0.5 
1.2 
1.4 


sys fixed to z spec = l.&9 
nearby part z phot ~ l.S 
nearby part z phot ~ l.S 


3.1 
3.2 

3.3 
4.1 

4.2 
4.3 


01 
04 
04 
04 
04 
04 


16:07.388 
16:08.461 
16:10.036 
16:07.398 
16:08.437 
16:10.051 


-24:04:01.62 
-24:04:15.53 
-24:04:32.56 
-24:04:02.01 
-24:04:15.53 
-24:04:33.08 


2.149 [2.130-2.323] 
2.324 [2.166-2.369] 
2.778 [2.759-2.814] 
2.199 [1.182-2.819] 
2.324 [2.166-2.369] 
2.244 [2.140-2.325] 


2.00 [1.61-3.74] 
1.99 [1.60-3.54] 


1.1 
0.5 
0.6 
0.4 

1 
0.5 


2.14 [1.59-3.97] 
1.81 [1.58-4.03] 


0.9 
2.3 

3 
2.3 
0.7 
0.9 


blended with 4.2 
blended with 3.2 


5.1 
5.2 
5.3 
C5.4 


04 
04 
04 
04 


16:07.773 
16:07.839 
16:08.043 
16:10.454 


-24:04:06.24 
-24:04:07.21 
-24:04:10.01 
-24:04:37.05 


2.362 [2.158-2.489] 
2.565 [2.469-2.733] 
2.166 [1.860-2.601] 


2.32 [1.74-3.63] 


1.4 
2.5 
0.6 
0.7 


2.54 [2.07-3.53] 


0.5 
2.6 
0.5 
0.6 


blended with member 
very plausible 


6.1 
6.2 
p6.2 
c6.3 


01 
04 
04 
04 


16:09.609 
16:09.946 
16:11.698 
16:11.422 


-24:03:42.64 
-24:03:45.31 
-24:04:21.12 
-24:04:14.95 


6.139 [0.804-6.662] 
1.470 [0.501-6.472] 


6.76 [3.81-7.31] 


2.0 
0.7 


6.29 [3.99-17.30] 


2.5 
0.4 


bimodal z pho( 
bimodal z phot 


c7.1 
c7.2 

p7.3 
c7.3 


04 
04 
04 
04 


16:09.552 
16:09.752 
16:11.166 
16:11.308 


-24:03:47.13 
-24:03:48.82 
-24:04:15.40 
-24:04:15.99 


2.376 [2.079-2.698] 
1.825 [1.801-1.982] 

2.378 [0.792-2.994] 


3.81 [3.36-4.63] 


3.2 
0.3 


4.18 [3.45-4.63] 


2.6 
0.7 




8.1 
8.2 
p8.3 


04 
04 
01 


16:08.783 
16:08.840 
16:11.266 


-24:03:58.05 
-24:03:58.83 
-24:04:24.70 


2.425 [2.266-2.622] 
2.414 [2.083-2.593] 


3.08 [1.32-3.41] 


0.5 
1.1 


2.97 [2.51-3.41] 


1.3 
0.3 




9.1 
9.2 
p9.3 
c9.3 


01 
04 
04 
01 


16:06.486 
16:06.605 
16:08.992 
16:09.149 


-24:04:42.90 
-24:04:44.78 
-24:05:00.00 
-24:05:01.23 


2.304 [1.922-2.547] 
2.386 [1.827-2.962] 


3.65 [2.87-4.33] 


0.8 
0.9 


3.93 [3.28-4.33] 


1.6 
0.2 


other cand. nearby 


10.1 
10.2 
pl0.3 
cl0.3 
cl0.3 


01 
04 
04 
04 
04 


16:06.244 
16:06.833 
16:09.158 
16:08.807 
16:09.818 


-24:04:37.76 
-24:04:47.12 
-24:05:00.45 
-24:05:01.94 
-24:04:58.69 


2.357 [2.236-2.526] 
2.356 [2.231-2.630] 
2.356 [2.231-2.630] 
2.124 [1.946-2.361] 
2.345 [2.115-2.775] 


4.63 [1.98-4.80] 


2.5 
0.9 


4.63 [3.04-4.80] 


1.4 
0.4 


mean chain z — 3.13 
other cand. nearby 


cll.l 
oil. 2 
cll.3 


04 
04 
04 


16:09.410 
16:09.196 
16:08.214 


-24:04:13.32 
-24:04:11.11 
-24:03:57.66 


0.717 [0.666-0.784] 
1.008 [0.833-1.092] 
0.619 [0.443-1.130] 


1.04 [0.99-1.19] 


2.7 

1.7 
1.2 


1.18 [1.17-1.20] 


1.8 
1.0 

1.2 




12.1 
12.2 
pl2.3 
C12.3 


04 
04 
04 
04 


16:09.230 
16:09.011 
16:06.966 
16:06.989 


-24:04:25.74 
-24:04:23.72 
-24:03:59.93 
-24:04:03.57 


1.618 [1.430-2.085] 
1.765 [1.031-2.043] 
1.288 [0.432-3.334] 


2.79 [1.94-3.41] 


0.6 
1.2 


1.82 [1.49-1.87] 


1.8 

1.8 




13.1 
13.2 
13.3 


01 
04 
04 


16:06.619 
16:07.711 
16:09.681 


-24:04:22.03 
-24:04:30.61 
-24:04:53.56 


3.255 [3.181-3.357] 
3.280 [3.008-3.481] 


4.35 [2.76-5.55] 


1.2 
1.2 
0.8 


3.24 [2.79-7.82] 


1.1 
3.7 
1.2 


mean chain z=3.83 


14.1 
14.2 
14.3 


04 
04 
04 


16:06.296 
16:07.450 
16:08.598 


-24:04:27.62 
-24:04:44.26 
-24:04:52.78 


1.765 [1.728-1.786] 
1.196 [1.093-1.243] 
1.773 [1.738-1.790] 


1.94 [1.60-2.42] 


2.8 
2.8 
1.2 


2.01 [1.81-2.27] 


1.1 

1.3 
0.9 


bright member nearby 


cl5.1 
cl5.2 
C15.3 


04 
04 
01 


16:06.292 
16:07.065 
16:08.560 


-24:04:33.67 
-24:04:42.90 
-24:04:55.38 


0.525 [0.205-1.774] 
2.905 [2.559-3.168] 


1.80 [1.39-1.95] 


2.1 
4.2 
3.0 


1.83 [1.64-1.95] 


1.5 
3.5 
1.5 


other cand. nearby 
other cand. nearby 


16.1 
16.2 
16.3 


04 
04 
04 


16:05.774 
16:06.799 
16:07.583 


-24 : 04:5 1 .22 
-24:05:04.35 
-24:05:08.77 


1 .964 [1.840—2 .079] 
1.803 [1.405-1.870] 


2.35 [1 . 77—3 . 16] 


0.8 
0.3 
0.5 


2.97 [2.66—3.41] 


1.9 
1.0 
0.6 




cl7.1 
cl7.2 
cl7.3 


01 
04 
04 


16:07.170 
16:06.866 
16:05.599 


-24:05:10.91 
-24:05:09.55 
-24:04:53.69 


2.467 [2.336-2.543] 


2.84 [1.92-3.81] 


0.9 
1. 1 
2.5 


3.87 [3.36-5.08] 


1.5 
1.9 

2.4 


part yields z phot ~ 3 


0I8.I 
cl8.2 
C18.3 


04 
04 
04 


16:06.258 
16:06.016 
16:07.416 


-24:05:03.24 
-24:05:00.06 
-24:05:12.28 


2.824 [1.314-3.083] 
2.766 [2.199-2.924] 


2.41 [1.61-4.47] 


0.4 
2.3 
0.6 


3.01 [2.87-4.72] 


0.8 
3.6 




C19.1 
cl9.2 
cl9.3 


04 
04 
04 


16:10.909 
16:10.777 
16:11.925 


-24:03:41.08 
-24:03:39.85 
-24:04:00.91 


2.774 [2.291-2.963] 
2.129 [1.229-2.611] 


3.01 [2.13-3.08] 


0.5 

1.8 
1.2 


3.60 [2.32-3.65] 


1.7 
3.8 
0.6 




c20.2 
c20.3 


04 
04 
04 


16:10.069 
16:10.478 
16:11.451 


' l n '4 III (; 'i 
-24:03:44.98 
-24:04:07.15 


1.306 [1.006-2.128] 
0.815 [0.469-4.155] 


2 56 [1 79—11 95] 


0.7 
2.6 




0.9 


bimodal z phot 


c21.1 
c21.2 
p21.3 
C21.3 


04 
04 
04 
04 


16:09.813 
16:09.865 
16:11.185 
16:11.047 


-24:03:46.67 
-24:03:47.32 
-24:04:11.82 
-24:04:07.73 


2.787 [1.868-3.039] 
2.446 [2.066-2.670] 

3.061 [2.771-3.244] 


4.47 [1.54-5.19] 


1.3 
0.3 


3.45 [3.20-5.19] 


1.1 
0.2 


moan chain 2 = 3.91 


C22.1 
c22.2 
c22.3 
p22.4 
p22.5 


01 
04 
04 
04 
04 


16:08.278 
16:08.204 
16:08.162 
16:08.508 
16:10.345 


-24:04:01.07 
-24:03:59.28 
-24:03:59.22 
-24:04:05.00 
-24:04:28.31 


0.362 [0.191-3.558] 
3.172 [0.511-3.427] 


1.74 [1.53-1.79] 


3.8 
2.3 
2.6 


1.67 [1.52-1.64] 


3.8 
2.3 
2.6 


cand. nearby 


C23.1 
c23.2 
c23.3 


01 
04 
01 


16:10.691 
16:09.505 
16:08.242 


-24:04:19.56 
-24:03:59.87 
-24:03:49.47 


2.410 [2.367-2.563] 
1.143 [1.109-1.214] 
2.057 [1.884-2.427] 


1.83 [1.23-2.45] 


1.6 
0.8 
1.1 


1.78 [1.16-2.54] 


1.0 
0.7 
1.8 


other cand. nearby 
other cand. nearby 



Note. — Column 1: arc ID. "c" stands for candidate and "p" for predicted location. For candidates the photo-z distribution, or 
identification by eye, was ambiguous; Columns 2 & 3: RA and DEC in J2000.0; Column J^: photometric redshift and 95% C.L.; Column 5: 
predicted and 95% CL. redshift by the eNFW model ; Column 6: reproduction distance of image from the observed location in the eNFW 
model; Column 7: predicted and 95% CL. redshift by the eGaussian model; Column 8: reproduction distance of image from the observed 
location in the eGaussian model; Column 9: comments. 



number depends also on the luminosity function and ob- 
servational depth, and on the mass shape parameters. 
However here our goal was to merely show that such a 
correlation exists and assess its order of magnitude. 

To test the effect of merger in M0416, we simulate two 
eNFW halos with similar parameters as in the resulting 
mass model including also the BCGs in their centers, ap- 
proaching each other on the line connecting them (Fig. 
Pfl). The observed separation between the two halos en- 



tails a ~ 20 — 60%(120)% increase in the total number 
of multiple-images(critical area), compared to the "far- 
away" initial position where each halo comprised its own 
(unmerged) critical curve. The exact amount is depen- 
dent on the true mass of each halo, including the true 
contribution of the galaxies (Fig. H. Note that as the 
two halos continue approaching eacn other, although the 
total number of multiple images may increase, the num- 
ber of multiple images per critical area, will decrease. 
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Fig. 3. — Number density of multiple images as a function of 
lens ellipticity (normalized to the circular case), for different con- 
centration and scale-radius parameters. In all probed cases, higher 
ellipticity boosts the lensing efficiency, generating more multiple 
images per critical area. We mark on the figure the measured 
effective ellipticity of M0416 and other (less) elongated clusters. 
Choosing those combinations of C200 and r s that yield comparable 
critical area to that of M0416, we obtain that the number density 
of multiple images observed in M0416 is ~ 2.5x higher due its 
elongation. 
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Fig. 4. — Critical area (Blue dash-dotted line) and number of 
multiple images versus the displacement between the two halos, 
normalized to a large-separation case. The Black curve shows the 
behavior when accounting only for the DM halos and BCGs, while 
the Red circle curve shows the same scenario but with more massive 
halos to account for the missing galaxies mass. In either case, it 
is clear that the merger contributes a boost of ~ 20 — 60% to the 
number of multiple images observed. 



To assess how extreme is the elongation of the criti- 
cal curves, we co mpare to Multidark/M USIC 2 numer- 
ical simulations (Sembolini et al. 2012). We used the 
80 clusters at z ~ 0.42 above the completeness mass 
limit M vir = 6 x lO 14 ft _1 M in a volume of l/i" 3 Gpc 3 . 



Each cluster was projected along 100 lines of sight. We 
searched for those projections where the critical lines 
both exceed 50" in at least one direction, and exhibit 
high elongations (axis ratio > 5). These conditions are 
satisfied in ~ 4% of the 8000 lens planes. 

4. SUMMARY 

We presented a SL study of M0416, in which we uncov- 
ered 70 multiple images and candidates of 23 background 
sources. We constructed mass models using two inde- 
pendent methods, both yield similar critical curves and 
mass profiles, in agreement also with independent larger- 
scale WL analysis. Compared to other well-known lens- 
ing clusters, M0416 exhibits a high number of multiple 
images for its critical area (~ 0.6 □'). 

We simulated the effects of ellipticity and merger on 
the lensing efficiency, and showed that (a) the number 
density of multiple images increases with ellipticity (the 
source-plane caustics get bigger, generating more multi- 
ple images for the same critical area), and (b) the critical 
area, and correspondingly total number of multiple im- 
ages, increase with lower separations between two merg- 
ing clumps, peaking when the two halos are either on top 
of each other or a few dozen arcseconds away, depending 
on their mass and shape. 

We conclude that the observed critical area size can be 
attributed to the merger, which boosts the mass in the 
center. For this given critical area, the high multiple- 
image number density can be explained by the observed 
elongation, which boosts the lensing efficiency by ~ 2.5 x . 
Backgroun d cosmic variance, estim ated typically at a ~ 
20% level rtSomeryille et al.||2004[ ) or double for z > 5 
(Trenti & Stiavclli 2008j), is likely to play only a small 
role in the increased number of multiple images, which 
are spread over a wide redshift range. 

This cluster shows once more that there exists a class 
of prominent lenses at redshifts around z ~ 0.4 — 0.5 
(and higher), probably due to their merging state and 
thus level of substructure and ellipticity which as we have 
shown, boost the SL properties. 
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